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ABSTRACT 

The 2-ni robotic Liverpool Telescope [LT) reacted promptly to the high-redshift (z = 4.048) gamma- 
ray burst GRB 060206. The afterglow was identified automatically and multicolor r'i' z' imaging 
program was triggered without human intervention. Combining our data with those obtained from 
later follow-ups provides a well-sampled optical light-curve from 5 minutes to >2d after the gamma 
event. The light-curve is highly structured with at least three bumps evident in the first 75 minutes, 
including a major rebrightening (Ar' « —1.6 at i « 3000s), interpreted as late energy injection. 
At early time {t ~ 440s), we find evidence for fast (Attest < 4s << f) variability, indicating on- 
going internal-engine activity. We emphasise that a low redshift GRB (z < 1) with similar intrinsic 
properties would have been interpreted completely differently due to undersampling of the light curve 
in the rest frame at early times; the light-curve behaviour of GRB 060206 should therefore not be 
considered peculiar. Finally, although the observed late-time steepening of the optical light curve 
resembles a jet break if taken in isolation, the lack of a corresponding change in the X-ray slope 
rules out a jet-break interpretation. Traditionally, GRB jet breaks have been inferred from optical 
data in the absence of simultaneous X-ray data. We suggest therefore that current estimates of the 
jet opening angle distribution might be biased by events like GRB 060206. Consequently, the GRB 
explosion energy distribution and event rates may have to be revised. 
Subject headings: gamma rays: bursts - cosmology: observations 



1. INTRODUCTION 



The Swift satellite (|Gehrels et all I2004D . launched 
in November 2004 and specifically designed to study 
Gamma Ray Bursts, is providing an overall rate of about 
100 real-time gamma localizations per year. A consid- 
erable fraction (w 75%) are reobserved within minutes 
by the on-board X-ray {XRT) and UV-optical {UVOT) 
narrow field-of-view instruments, whilst an increasing 
number of ground-based robotic facilities are participat- 
ing actively in the worldwide follow-up program by re- 
sponding rapidly to satellite alerts. As a result, the 
number of GRBs with early-time light curves and spec- 
troscopically determined redshift has increased dramat- 
ically over the last 18 months. Thus, in addition to 
important recent results such as the detection of short 
burst afterglows and the observati on of the most dis- 
tant (z>6) GRB obse rved to date ijHaislip et alJl2005t 
lAntonelh at al. ll2005(l ). the following facts are emerg- 
ing: aj a higher meaii redshift {z « 3) for Swift GRBs 
l|jakobsso n et aljl2006|) c ompared with thos e detected in 
the pre-Swift era (z sa 1 ijBerger et al.ll2005j) ): h) 990123- 
like (^kcrlof et al. 199^ prompt flashes are on average 
weaker and less frequent than expected; c) & popula- 
tion of faint bursts (truly dark only in absence of an 
early, deep follow-up) providing usefu l constraints on ex- 
isting opti cal suppression models ( Homi ng et alJ 120051: 
lOates et"al.. 20061: d) huge, rapidly rising and decaying 
flares observed in the X-ray light-curves, indicating prob- 
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able late-time activity of the inner engine l|Burrows et alJ 
^005) . The growing evidence for continued activity of the 
central engine, the detailed physics of which remains elu- 
sive, and the opportunity to probe the interaction of ex- 
panding shells with the circumburst/interstellar medium, 
particularly at high redshift, demand dedicated prompt 
foUowup observing programs. Large robotic telescopes 
such as the 2-m Liverpool Telescope ar e optimised for 
such foUowup and are proving effective l)Guidorzi et alJ 
^05, 2006). 

Early afterglow light curves exhibit a considerable va- 
riety and rebrightcnings (hereafter called humps) are 
becoming more frequently observed. GRB 021004 

jiBersier e t 311120031 iMatheson erani2003t .Ho lland et alJ 
2003) remains one of the best studied events [z = 
2.335) and a number of different models have been 
invoked to explain its temporal beha viour: vari- 
able density profile fLazzati et al. 2002), passage of 
the break freauen cv (.Kobavashi fc Zhang 20 03. ), re- 
freshed sho cks (Na kar et alJl2003HBiornsson et al.ll2004l: 

Ide Ugart~Postigo et al.ll2005t > . angular dependence of 
the energy profile on the jet structure ( Nakar et al. 2003) . 
Another interesting case is GRB 030329 [z = 0.17), 
preferably modeled according to t he refreshed sh ocks 
(late energy-injections) scenario ifHuang et al.ll200^ . In 
this case, initially slower shells catch up with the exter- 
nal shock front, which, in the meantime, is decelerat- 
ing due to the interaction with the surrounding medium. 
The early afterglow of GRB 050502a {z = 3.793) ex- 
S£fiik}il5'4i*jSffle3th decay, well- fitted by a power-law with 
mdex a = 1.2 (F^, oc t~"), followed by a an achromatic 
bu mp emerging at tre.st ~ 600s that has been interpreted 
bv iGuidorzi et all l)2005j) as due to density clumps in a 
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uniform medium. More recently, an increasing number 
of events, in p artic ular at high-z, exhibit rebrightenings 
(jRvkoff et al.il200l Istmek et al.ll2006|) . These observa- 
tions support the hypothesis such behaviour might not 
be unusual. We report here multicolor robotic observa- 
tions of GRB 060206 carried out by the Liverpool Tele- 
scope (LT) and revealing a structured behaviour of the 
early lightcurve. The 2-m aperture of the LT, the opti- 
mal observing con ditions and a pre cise SDSS pre-burst 
field calibration l|Cool et al.l l2006(l contributed to the 
remarkable quality of t he dataset. Timely later pho- 
tometry irom RAPT OR l|Wozniak et al.l2006|) and MDM 
ijStanek et al.ll2006(l allowed us to compile an exception- 
ally information-rich light-curve catching the different 
phases of the early afterglow. 

2. OBSERVATIONS AND ANALYSIS 

On 2006 Feb 06, 04:46:53 UT^ Swift-BAT detected 
GRB 060206 (trigger 180455) at Galactic Coordinates 
l=78.07deg, b=78.28deg4 (3 arc-min 90% containment 
radius). In the 15 — 350keV band the prompt gamma 
profile is well described by a single quasi-gaussian peak 
with duration 7s±2s. Swift-XRThegan observing 65s af- 
ter the BAT trigger time a nd entered the South Atlantic 
Anomaly soon afterwards ijMorris et al.ll2006|l . UVOT 
started the on-target monitoring practically at the same 
time. As usual, a number of ground-based facilities were 
activated to follow-up the event. The redshift ha s been 
spectroscopically measured to be 2; = 4.048 ( Fvnb o et al.l 
I2QO6). Robotically triggered photometric observations 
with the Liverpool Telescope began a.t t = 309s af- 
ter the GRB, at an airmass of 1.01. A sequence of 
2 X 3 X 10s r' exposures, named detection mode, aims 
at finding potential afterglow candidates and automati- 
cally triggers the best suitable early follow-up observing 
pr ogram. Detai ls of the procedure are reported elsewhere 
bv lGuidorzi"er al. ( 2006). In the GRB 060206 case, a re- 
peated sequence of r'i'z' 120s exposures was initiated. 
The afterglow is detected with high S /N in all the single 
exposures. Calibration was facilitated and greatly im- 
proved by the ava ilability of SDSS pre-burst photometry 
ijCool et al.l2006j) . We adopted 5 calibrated star-like PSF 
field objects to adjust the zero-point of the single images. 
The zero points were stable throughout the observation 
sequence, and consistent with nominal LT values^. Pho- 
tometry was carried out independently using the Star- 
link/GAIA and DoPHOT programs. The datasets are 
in complete agreement with each other, and the results 
have been combined to further reduce the systematics. 
Magnitud es are converted into flux densities F^, (mJy) 
foUowing iF^ik^igita et alJ l)1996l ^. Results are summarized 
in Tabled 

Figure^shows the complete lightcurve of GRB 060206. 
100 RAPTOR (Wozniak ct al. 200^ and 101 MDM 
ifStanek et alJ ^2006) photometric points have been in- 
cluded together with selected GCNs. The RAPTOR 
unfiltered values, calibrated with respect to USNOBl.O- 

^ t=0 througho ut the paper. 

4 According to ISchlegel e rai? fToo l) E(B -V ) = 0.013. The 
estimated extinction, from ICardclH c t alJ 1198911 curve, is A^/ = 
0.036, Ai, = 0.027 and A^, = 0.019. 

^ Confirmed also by the GMT in La Palma. 
Nightly av eraged atmospheric extin ction: r'=0 .076. 
http: //www. ast.cam.ac.uk/ ~dwe/SRF/camc_extinction.html| 



R2 magnitudes ijWozniak et al.l r2006j) . have been flux- 
recalibrated with respect to the earlier LT data. This is 
equivalent of shifting the RAPTOR data by -1-0.3 magni- 
tudes. This approach relies on a relatively stable colour. 
At a first glance this assumption might look uncertain 
since we are actually presenting evidence of a spectral 
evolution across the bump. However, by assuming a re- 
alistic A{r' — i')tump ~ 0.4, a rather small Ar'^y^ « —0.07 
is estimated. This systematic correction has been added 
to the synthetic r' late photometry data. In any case this 
isn't in any way affecting our c onclusions. The MDM 
data, as in iStanek et alJ l|2006|) . have been shifted by 
4-0.22 magnitudes with respect to RAPTOR and accord- 
ingly converted into physical units. XRT dsda have been 
processed with the ftool 'xrtpipeline' to produce cleaned 
level II event files, applying default screening constraints. 
We select only events with grade for both PC (Pho- 
ton Counting) and WT (Windowed Timing) modes. For 
PC mode, we initially considered a 20 pixel-radius cir- 
cle centred on the source. From the PSF profile, it was 
clear the data were pa rtially affected by pile-up; following 
iVaughan et al.l l)2006|) . we corrected for this by exclud- 
ing photons within a 5-pixel radius inner circle and then 
renormalised the resulting light curve for the fraction of 
the PSF profile considered. Background spectrum and 
light curve have been taken from four 50-pixel radius cir- 
cles with no sources. No significant pile-up was found 
to affect WT data, fo llowing the procedure described in 
iRomano et al.l l)2006j) . We adopted empirical ancillary 
files for the spectra created with the ftool 'xrtmkarf. 
We used the latest spectral response matrices in the Cal- 
ibration Database (CALDB 2.3). 

2.1. r' Lightcurve Fit 

The analytical model adopted for the fit, up to the 
break at ~ 3 x 10^ s, is: 

Fjt) = y F, • (1) 

This is the sum of a number of smooth ly connected 
broken power-laws ijBeuermann et anil999|) . ai and a2 
depend on j and are the exponents a,tt << tj and t >> tj 
and tj the bump (break) time. When ai < and 012 > 0, 
a bump is produced. A typical light-curve break, on the 
other hand, is obtained when ai,a2 > and 012 7^ ai. 
The constant n controls the smoothness of the transi- 
tions, but has a relatively small influence in this context 
(see Table [JJ . Fj is the individual flux contribution at 
the bump (break) time. In the case of GRB 060206 the 
j index runs from to 5 (see Figure ^1. The j = 
segment is, by definition, associated with tj = and 
is therefore a simple power-law fitting the earliest LT 
photometric points. The j = 1 segment accounts for 
the sharp, earliest bump occurring at t « 440s and the 
subsequent decay. The next segment describes three r' 
points lying at t « 2000s. However, this fit at 2000s 
is not unique as one may also fit the flattening of the 
curve around this time with Aa w —0.8, beginning im- 
mediately after t = 1200s. Although this scenario may 
seem simpler than introducing a bump component via 
J = 2, assuming an achromatic flattening in the r' band 
curve for t — 1200s : 1900s (constrained by r' sampling) 
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would in turn imply a rebrightening in z' and i' . This 
confirms that a bump is required in this time interval. In 
our analytical model the bump time is arbitrarily fixed 
to t2 — 1500s, and the shape fixed by the earliest case 
template. The only variable is the peak flux Fi that is in 
any case non significantly affected by ti- In other words, 
a small rebrightening at t w 1500s is necessary because it 
is not possible to fit the light curve in all 3 colors {r'i'z') 
without including a bump at that time. The next seg- 
ments, j = 3,4: provide a good fit to the biggest double- 
peaked bump (<3 = 3529s ± 26s, = 5234s ± 79s). Most 
importantly, the cosmological time dilation (l+z) implies 
that the intense activity described so far took place in less 
than 15 minutes in the rest frame. The last feature be- 
fore the final break, possibly associated to a density vari- 
ation, takes place at t w 16000s and is well described by 
the j — 5 term. In absence of a systematic photometric 
sampling close to the late break {t ~ 0.6d), we fitted the 
t > Id MDM points with a simple, independent power- 
law. The analytical model validity is thus limited to 
t < 50000s, before the occurrence of the inferred break. 
The result is a steepening from ai = 0.95 {t « ttreak) 
to a2 « 1.8 {t >> tbreak)- Wc notc that these values 
differ slightly from those found by[Stanck ct al. (2006) 
who derived ai ~ 0.7 to a2 ~ 2.0. The difference is to 
be attributed probably to the j — 5 bump that is recog- 
nized as a separate feature i n our view, but inclu ded in 
the general pre-bump slope in lStanek et al.l l)2006j) . That 
late break nature is discussed in more detail in the next 
Section. 

The analytical model fit was performed separately in 
three different time ranges: 300 : 2500s (6 Degrees Of 
Freedom), 2500s : 15000s (126 DOF) and 7000s : 50000s 
(74 DOF). Since the fit parameters are common, the 
procedure is iterated to let them converge to the same 
values in the three intervals. Convergence is achieved 
after ^ 10 iterations, and the result is then stable in- 
dicating that a true global minimum has been reached. 
A summary of the parameters resulting from the fit is 
reported in Table [3 In Figure the single components 
are shown together with the resultant fitted curve. 

2.2. Multi-epoch Spectral Energy Distribution Analysis 

Figure 121 shows the pre-bump rest-frame SED interpo- 
lated at < = 785s, where the r' lightcurve is relatively 
smooth, from LT and XRT data. The XRT spectral 
shape at this epoch have been assumed coincident with 
the PC early template {t = 141s : 341s). The X-ray flux 
att = 785s, is derived from an a = 1 (Tabled power-law 
extrapolation. This seems a quite reasonable assumption 
(see Figure Further, no significant differences are 
seen between the XRT spectra extracted from PC and 
WT modes. An extrapolated UVOT-B point has been 
added to provide important qualitative constraints. The 
average time behaviour of the B light-curve is assumed to 
follow the well-sampled r' case. An Fr' {t) integral is per- 
formed over the UVOT time span (t = 1662s : 34946s, 
iBovd et al.l l|2006(l ). The averaged F^r'' is directly com- 
pared to Fg^'-^^ . The following simple ratio relation pro- 
vides an estimate of Fb corresponding to the pre-bump 
epoch of interest: 

F^r'' : Fr' (T) = F^^O"^ : Fb{T) (2) 



At t = 785s the four optical points are not reconcilable 
with a single power-law. The r' flux is slightly affected, 
according to iFvnbo et alJ l)2006f) . by the deep Damped 
Lyman Alpha DLA absorption feature. The i' and z' 
values, on the other hand, can be very well extrapolated 
up to I'rest ~ ASAOPHz following a rather typical /3 = 0.8 
slope (F^ oc i/-'3). The UVOT - B point qualitatively 
confirms the pronounced HI absorption and is not in- 
cluded in the fitting procedure. Even if sparsely sam- 
pled, the positions of the i' and z' photometric points 
in the time interval t — 1000s : 3000s clearly indicate ei- 
ther a SED evolution or an additional unresolved activity 
(or both combined). A rigid shift^ of the r' light-curve, 
interpreted either with the bump or the flattening, indi- 
cate that the i' and z' fluxes are significantly higher than 
expected. According to the photometric errors, the dis- 
agreement significance has been estimated in 4(7 in the 
bump scenario (i', 9cr for z'). Analogous results assum- 
ing an achromatic flattening (5cr and llcr). We suggest 
that the SED evolution scenario is more conservative, 
but caution that a generally accepted quantitative pic- 
ture would require better sampling of the light curves. 
As an additional note, we remark that the simple rescal- 
ing described in footnote is again able to account for the 
z' points calculated during the big bump rising slope. 
This might be a possible indication, even if not conclu- 
sive, that the SED is already relaxed before the bump 
reaches its peak. 

In Figure El we show the post-bump SE D at t = 5680s 
based on the FAIRITIEL JHK detections CAlatalo et all 
2006) {tstart = 2.91/i) and on the earliest MDM R point"^. 
Both are back-extrapolated according to an a = 0.95 
temporal slope (table (SJ. At this late time, already 
on the smooth side of the afterglow, the IR-X fluxes 
{vrest = 0.7 : 8000Pi/z) are well fitted by a common 
13 = 0.93 ± 0.02 power-law [x^/d.o.f. = 3.46/10). How- 
ever, a broken power-law with fixed indices (3opt = 0.7, 
Px = 1-2 is not ruled out {x^/d.o.f. = 10.04/10). The 
latter is considered in the light of the interpretative 
model discussed later. 

Post-break IR photometry l|Terada et al.l I2006D and 

MDM late observations give Pir^o = 0.7 ± 0.3. A sum- 
mary of the spectral results is given in Table O 

3. DISCUSSION 

3.1. Re- Brightening 

After the remarkable rebrightening at ^ Ihr and be- 
fore the break around t ~ 5 x lO'' s, the optical flux 
is described by a single power law with a ~ 1, except 
small wiggles. This indicates that the blast wave radi- 
ating the afterglow had a significant transition from one 
Blandford-McKee (BM; Blandford&McKee 1976) solu- 
tion to another BM solution around t — Ihx. Since blast 
wave solutions depend only on two parameters (explo- 
sion energy and ambient density) possible scenarios are 
energy injection (Rees& Meszaros 1998; Kumar& Piran 
2000; Sari & Meszaros 2000) and a density-jump medium 

^ The multiplicative factors for i' and z' are determined ac- 
cording to the earliest epoch and are respectively 1.50 ± 0.05 and 
1.72 ±0.05. 

^ We adopted the original lStanek et alJ I200fii) R flux instead of 
the estimated r'. The SED is thus completely unaffected by the 
LT-to-RAPTOR conversion systematics described in detail in the 
text. 
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(Ramirez-Ruiz et al. 2001; Dai & Lu 2002). The degen- 
eracy is broken by the X-ray observations. Since there 
is no signature of the coohng break in the decay phase 
of the X-ray afterglow (Morris et 2006), the X-ray band 
is likely to be located above the cooling frequency at 
t ~ 1 hr. This implies a high ambie nt density, consistent 
with spectroscopic determinations l)Fvnbo et al. 2006). 
The late time break in the optical light curve also favors 
'^x > at ~ 1 hr as we will discuss below. If the X-ray 
band is above the cooling freque ncy, the flux does not 
depen d on the ambient density ijFreedman fc WaxmanI 
120011) . We can conclude that the X-ray rebrightening is 
not due to a jump in the ambient density. 

According to Kumar & Piran (2000), in the absence of 
a density enhancement, the flux depends on the explosion 
energy as: 



(3) 
(4) 



Since the optical flux increased by a factor of ~ 4, the 
rebrightening indicates there was a huge impulsive en- 
ergy injection AE ^ 1.8i?o at ~ 1 hr where Eq is the 
blast wave energy right before the rebrightening. A typ- 
ical p=2.3 value is assumed for the power-law distribu- 
tion of electrons energies. Consistently with observations 
(Figure^l , the X-ray flux is expected to increase by a fac- 
tor of 3. A rather smooth transition from one solution to 
the other at most wavelengths is expected. However, this 
is not the case at frequencies for which there is significant 
emission from the reverse shock (Kumar & Piran 2000; 
Kobayashi & Sari 2000). The sharp rise of the rebright- 
ening implies that the typical frequency of the reverse 
shock is close to the optical band. The small wiggles 
before or after the major rebrightening can be explained 
by additional minor energy injection or density clumps 
in the ambient medium. 

The energy injection is probably due to a decrease in 
the Lorentz factor of the outflow toward the end of the 
prompt GRB. This slow outflow collides with a blast 
wave when the blast wave have slowed down as a result 
of sweeping up the ambient material. Energy injection 
could also be caused by a long lasting activity of the 
central engine (Burrows et al. 2005). The energy in the 
blast wave at late times (> Ih) is larger than that at the 
deceleration time by more than a factor of ~ 2.8. This 
requires the prompt gamma-ray emission process (inter- 
nal shock s) to be more efficient than previously consid- 
ered (lokaet al. ll2005HNousek et 31.1.2006; .Zhang et all 
l20f)6tlGranot. Konid fc Piran II2006D . 

3.2. Jet Break? 

The optical light curve of GRB 060206 steepened from 
a = 0.95 ±0.02 to 1.79 ±0.11 around t = (3 ~ 9) x 10^ s. 
Although this steepening resembles a jet break, the X- 
ray light curve from t = 5 x 10"^ s to 7 x 10^ s re- 
mains consistent with a single unbroken power law de- 
cay with a = 1.35 ± 0.15 (Morris et al 2006). This 
clearly contradicts the monochromatic break predica- 
tion by the jet mode l. In the wind model p{R) oc 
ijChevalier fc Lll999() . the cooling frequency increases 
as a function of time Vc t(3s-4)/(8-2s)_ rj,]^^ passage 
of the cooling frequency can produce a steepening of 



5a = (3s — 4)/ (16 — 4s) in the optical band first and 
in the X-ray later. The observed steepening 5a — 0.84 
gives s ^ 2.7. The cooling frequency rapidly increases 
as t^'^, and it will cross the X-ray band around ^ 10^ s. 
The observed unbroken power law X-ray light curve up to 
t ^ 7 X 10^ s is consistent. However, this scenario does 
not account for the difference of the decay indices be- 
tween the pre-break optical (ttopt = 0.95) and the X-ray 
{ax = 1.35) light curves, which should be on the same 
spectral segment and the decay rates should be identical 
to be consistent with the scenario. 

A possible solution is that a blast wave initially prop- 
agates into a constant density medium p = rupn, and 
around i ~ 3 x 10"* s it breaks out in a wind-type p oc 
medium. Using the isotropic gamma-ray energy about 
6 X lO^'^ergs (Palmer et al. 2006), the transition radius 
is estimated in equation (5) in which ( = (I + z)/5 and 
ni — n/10 protons cm"^: 



Rtricm) - 3xl0"C"^^^(: 



t 



E 



1/4, 



'3xl04s' ^QxlQ^'^erg' 

(5) 

Although a large amount of energy is injected to the 
fireball after the prompt emission, the energy depen- 
dence of Rtr is rather weak. In this scenario, the or- 
der of the break and observational frequencies remains 
as Urn < Vopt < Vc < Vx after the major re-brightening 
t > \ hr. During the constant medium phase, the time 
and frequency dependences of the optical and X-ray flux 
are given by Equations (6) and (7) respectively (Sari, 
Piran fc Narayan 1998): 



F.OC i-(3f-2)/4^-p/2 



(6) 

(7) 



With the standard value of p = 2.3 these give rea- 
sonable fits to the observational decay indices, while the 
estimated optical spectral index (3 ~ 0.7 is slightly shal- 
lower than the observations j3 ^ 0.93 ± 0.02. Moreover, 
a spectral break is not directly inferred between optical 
and X-ray that seem to lie on a single slope. However, as 
explained in the previous Section, the predicted (3 ~ 0.7 
value is not significantly ruled out when a broken power- 
law is invoked. When the blast wave propagates into 
a wind medium, as the ambient density decreases, the 
optical flux drops as Fopt oc i^"™ where 



[2s + 3(p-l)(4-s)] 
(16 -4s) 



1.6 



(8) 



for p=2.3 and s=2.5. Since the X-ray band is above 
the cooling frequency, its flux does not depend on the 
ambient density, and the X-ray decay index does not 
change before and after the break-out. Considering the 
scaling between the observer time and blast wave radius 
t oc i?^''"'*^ , the steep segment of the optical light curve 
requires that the ambient density should behave as R~^-^ 
between Rtr and iRtr- The steepening from a = 0.95 
to 1.79 in the optical light curve looks like a typical jet 
break. The X-ray light-curve, however, shows no sign of 
a change in decay slope. This conflicts with the jet break 
potential interpretation based on the optical light-curve. 
In most previous events. X-ray data are not available at 
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the "jet breaks", and the breaks are interpreted based 
on optical data. The current estimate on the jet opening 
angle distribution might be biased by events like GRB 
060206. Consequently, the GRB explosion energy distri- 
bution and the event rates may need to be reviewed. 

4. CONCLUSIONS 

We presented a comprehensive study of the high red- 
shift GRB 060206 based on early robotic observations, 
later published photometry, GCN circulars and Swift 
public data. The composite monochromatic light-curve 
is exceptionally well sampled and extremely structured. 
An early, A< << t rebrightening suggests internal en- 
gine activity (internal shocks) continuing well beyond the 
gamma prompt emission. Later on, a high-significance 
SED evolution is inferred from multi-color r'i' z' pho- 
tometry carried out by the Liverpool Telescope; this is 
probably associated with another pre-bump rebrighten- 
ing. One hour after the GRB a significant, steep double- 
peaked rebrightening occurred, increasing the flux by a 
factor ~ 4 with respect to the pre-bump values. We in- 
terpret the cause of this rebrightening to be due to a 
late-time energy injection. The energy in the blast wave 
at late time {t > Ih) is a factor of 2.8 larger than that at 
the deceleration time, implying a higher internal shock 
efficiency than previously considered. Following the big 
bump, a rather typical spectral/temporal afterglow be- 
haviour is observed, with a possible density variation fea- 



ture superimposed on the smooth power-law decay. We 
interpret the fact that the steepening in the optical light- 
curve observed at t ~ 0.6d post-burst is not matched by 
similar behaviour in the X-ray as due to a transition in 
the density distribution of the ambient medium, likely 
from a constant density medium to a wind medium. In 
the absence of X-ray data the observed optical steepen- 
ing would have been interpreted as a typical jet-break. 
Since most of the so far observed jet-breaks lack cor- 
responding X-ray confirmation, we suggest that current 
estimate on the jet opening angle distribution might be 
biased by events like GRB 060206. The GRB explosion 
energy distribution and the event rates are thus poten- 
tially affected. 
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Monfardini et al. 




Fig. 1. — Complete GRB 060206 light-cur ve. Ea rly (t < 3600^) mu lt i-color r'i'z' ph otomet ry from Liverpo ol Tele scope robotic 
observ ations. Later data from IWozniak et alj 120061) . IStanek et all 120061) . IMalesani etal] 120061) . jReichart eTall II2006I) . IHaislip etaT] 

120061) . Also shown: analytical model best fit, separate contributions (dotted lines), Jfi?T counts rates (right-side axis). 



GRB 060206 Early-Time Light Curves: GRB Energetics and Environment 



1000 
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10 



UVOT« 
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Post-bump X 0.1 't--*-- 



0.01 I ^ • • ^ • • ^ • • ^ • • 
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Rest Frame Frequency (PHz) 

Fig. 2. — GRB 060206 pre- and post -bump Rest-frame SED {t = 785s and t = 5680s). Optical points from Liverpool Telescope robotic 
observations, IS tanek et all 120061) and IBovd etaT] 1200611 . The 5680s fluxes are divided by 10. Dotted lines: = 0.80 and 13 = 0.93 
power-law fitting respectively pre- and post-bump points. In the text we report evidence of a significant SED evolution occurring in 
t = 1000s ; 3000s. The fit results are summarized in Table ISl 
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TABLE 1 

Liverpool Telescope GRB 060206 Optical Photometry 







rijxposure 


iVlagiiitucle 


r iux 




Col 








r' 


319 


10 


16.93 ±0.02 


0.614 ±0.011 


r' 


341 


10 


17.02 ±0.02 


0.565 ± 0.010 


r' 


363 


10 


17.04 ± 0.03 


0.555 ± 0.015 


r' 


415 


10 


17.12 ±0.03 


0.515 ± 0.014 


r' 


437 


10 


16.95 ±0.02 


0.603 ± 0.011 


r' 


459 


10 


16.96 ±0.04 


0.597 ±0.021 


r' 


574 


30 


17.12 ±0.02 


0.515 ± 0.009 


i' 


785 


120 


17.04 ± 0.04 


0.555 ± 0.020 


z' 


885 


120 


17.03 ±0.03 


0.560 ±0.015 


r' 


1036 


120 


17.80 ±0.02 


0.275 ± 0.005 


r' 


1169 


120 


17.93 ±0.03 


0.244 ± 0.007 


i' 


1336 


120 


17.23 ± 0.03 


0.466 ± 0.013 


i' 


1487 


120 


17.24 ±0.05 


0.461 ± 0.021 


z' 


1762 


120 


17.06 ± 0.03 


0.545 ± 0.015 


r' 


1921 


120 


17.98 ±0.03 


0.233 ± 0.007 


r' 


2052 


120 


18.02 ±0.03 


0.225 ± 0.006 


r' 


2185 


120 


18.02 ±0.03 


0.225 ± 0.006 


i' 


2355 


120 


17.31 ± 0.04 


0.433 ± 0.016 


i' 


2487 


120 


17.29 ± 0.04 


0.441 ±0.016 


i' 


2619 


120 


17.31 ± 0.03 


0.433 ± 0.012 


z' 


2807 


120 


17.02 ± 0.04 


0.565 ± 0.020 


z' 


2939 


120 


16.85 ± 0.04 


0.661 ± 0.024 


z' 


3071 


120 


16.41 ± 0.04 


0.991 ± 0.036 


r' 


3273 


120 


17.04 ± 0.04 


0.555 ±0.020 


r' 


3405 


120 


16.82 ± 0.02 


0.679 ± 0.012 


r' 


3538 


120 


16.73 ± 0.02 


0.738 ± 0.013 


r' 


3670 


120 


16.70 ± 0.02 


0.759 ± 0.014 



^Timc delay with respect to the GRB trigger time, to = 
0.09302 UT. 
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TABLE 2 

GRB 060206 SDSS-R lightcurve fit results'''''^ 



Segment 






0-2 




Comment 




(mjy) 






(s) 







n.d. 


n.d. 


0.94 ±0.17 





a2=early slope 


1 


0.160 ±0.016 


-30 (fixed) 


1.2 ±0.5 


439 ±3 




2 


0.076 ± 0.007 


-30 (fixed) 


1.0 (fixed) 


1500 (fixed) 




3 


0.566 ±0.020 


-8.0 ±0.5 


4.3 ±0.6 


3529 ± 26 


big bump 


4 


0.356 ±0.027 


-4.9 ±0.7 


0.95 ±0.02 


5234 ± 79 


<^2=post-bump slope 


5 


0.029 ± 0.002 


-3.7±0.7 


3.4 ±0.5 


16037 ± 397 


density bump 


6 


n.d. 


n.d. 


1.79 ±0.11 


Si 53000 


012= slope after final break 



^Refer to the text following equation 1 for an explanation. 

'^The n exponent is fixed to 2.5 for all the segments. 

"^The segment 6 has been fitted separately with a simple power-law. 



TABLE 3 

Fitted spectral slopes from multi-epoch SED (Figure and text 



Timers) 






Po-x 


Comments 


785 


1.42 ±0.58'' 


1.15 ± 0.55*^ 


0.81 ± 0.03 


X?>_xM°-/ = 4-26/4 


RJ 1800 


2.05<= 


1.15 ± 0.55*^ 


n.d. 


PoT^Px 


5680 


0.84 ±0.14^ 


0.84 ±0.15 


0.93 ± 0.02 


X?>_xMo-/ = 3-46/10 


Si 108000 


0.7±0.3<= 


n.c. 


n.c. 


X?fl_oMo-/ = 0-85/1 



^Excluding B (UVOT) value affected by HI absorption 
^Assumed fromXi?T spectrum in t — 141s — 341s. See text 

^^Calculated from r' i' z' points interpolated according the additional bump hypothesis 
(model-dependent) 
^From back-extrapolations in JHK (PAIRITIEL) and R (MDM) 
'^SUBARU-IR & MDM-R 



